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Abstract. We propose and develop the Bethe Ansatz method for the Heun equa- 
tion. As an apphcation, holomorphy of the perturbation for the BCi Inozcmtsev 
model from the trigonometric model is proved. 



1. Introduction 

Olshanetsky and Perelomov proposed a family of integrable quantum systems, 
which is called the Calogero-Moser-Sutherland system or the Olshanetsky- Perelomov 
system [^. In the early 90's, Ochiai, Oshima and Sekiguchi classified integrable mod- 
els of quantum mechanics which are invariant under the action of a Weyl group with 
certain assumptions [Q. For the case of {N > 3), the generic model coincides 
with the BCn Inozemtsev model. Oshima and Sekiguchi found that the eigenvalue 
problem for the Hamiltonian of the BCi (one particle) Inozemtsev model is trans- 
formed to the Heun equation with full parameters 0. Here, the Heun equation is a 
second-order Fuchsian differential equation with four regular singular points. 

In this paper we will investigate solutions to the Heun equation motivated by the 
analysis of the Calogero-Moser-Sutherland systems. 

In the papers holomorphy of the perturbation for the Calogero-Moser- 

Sutherland model of type Ajy from the trigonometric model was proved, and in |^ 
some results were obtained by applying the Bethe Ansatz method. Note that the 
Bethe Ansatz for the model of type was established by Felder and Varchenko |^ 
by investigating asymptotic behavior of an integral expression of a solution to the 
KZB equation. 

For the case of type Ai, the Bethe Ansatz method was found in the textbook of 



Whittaker and Watson |T^, although it was covered in the chapter entitled "Lame's 
equation" , and the phrase "Bethe Ansatz" was not used. Note that the Lame equation 
is a special case of the Heun equation p. 

In this paper we will propose and develop the Bethe Ansatz method for the Heun 
equation, and holomorphy of the perturbation for the BCi Inozemtsev model from 
the trigonometric model will be obtained. 
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The method "Bethe Ansatz" appears frequently in physics. In particular the Bethe 
Ansatz method has great merits for investigating solvable lattice models |^. In this 
paper, we use "Bethe Ansatz" in a somewhat restricted sense, similar to the usage 
in 0, 1^. The Bethe Ansatz method replaces the problem of finding eigenstates and 
eigenvalues of the Hamiltonian (for example ( |2.2D ) with a problem of solving transcen- 
dental equations for a finite number of variables which are called the Bethe Ansatz 
equations (for example ( p.29D ). In our case, we use elliptic functions to describe the 
Bethe Ansatz equation. 

It would be very difficult to solve the transcendental equations explicitly. Instead, 
by applying the trigonometric limit the transcendental equations tend to algebraic 
equations, which may be more easily solvable. 

On the other hand, the Hamiltonian of the BCi Calogero-Sutherland model appears 
after applying the trigonometric limit, and the eigenstates for the BCi Calogero- 
Sutherland model are described by Jacobi polynomials. Thus the Jacobi polynomials 
and the solutions to the Bethe Ansatz equation for the trigonometric case are closely 
connected, and they are related to the perturbation. In fact holomorphy of the 
perturbation for the BCi Inozemtsev model from the trigonometric model follows 
from holomorphy of the solutions to the Bethe Ansatz equation with respect to a 
parameter of a period of elliptic functions. 

Now we examine the perturbation. Based on the Hamiltonian Hqs of the BCi 
Calogero-Sutherland model, the Hamiltonian H of the BCi Inozemtsev model can be 
regarded as a perturbed one, i.e. H = Hcs + Vo + Y^°l^ Vi{x)p\ where Vq is a constant 
and Vi{x) {i G Z>i) are functions. We can calculate eigenvalues and eigenstates of 
H as formal power series in p by the perturbation. For our case, holomorphy of the 
perturbation can be shown and it implies convergence of the formal power series in p 
if \p\ is sufficiently small. 

This paper is organized as follows. In section H, we clarify the relationship between 
the Heun equation and the BCi Inozemtsev system [0 . In section |^, we introduce 
the Bethe Ansatz method for the BCi Inozemtsev system, or equivalently for the 



Heun equation. In section 3.1, invariant subspaces of doubly periodic functions are 



introduced. In section |3]^, we obtain an integral representation of eigenfunctions of 
the BCi Inozemtsev system. In section |3.3| , we introduce the Bethe Ansatz method for 
the BCi Inozemtsev system and describe the Bethe Ansatz equation. In section 
we rewrite the Bethe Ansatz equation in terms of theta functions in order to consider 
the trigonometric limit. In section ^, we consider the trigonometric limit and solve the 
trigonometric Bethe Ansatz equation. In section ^, relationships between the Bethe 
Ansatz method and the spectral problem on L^-space for the BCi Inozemtsev system 
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are clarified. Then holomorphy of tlie perturbation for tlie BCi Inozemtsev model 
from the trigonometric model is obtained. A key observation is that holomorphy 
of the perturbation follows from holomorphy of the solutions to the Bethe Ansatz 
equation. In section ^, we give comments, and in the appendix, we note definitions 
and formulae of elliptic functions. 



2. The Heun equation and the Hamiltonian of the BCi Inozemtsev 

SYSTEM 

The expression of the Heun equation in terms of elliptic functions was already 
known by Darboux in the 19th century, and this appeared in connection with the 



soliton theory by Verdier and Treibich [|TT|. In the relationship between the 
eigenfunctions for the BCi Inozemtsev model and the solutions to the Heun equa- 
tion was demonstrated. In this section, we recall the relationship between the BCi 
Inozemtsev model and the Heun equation. 

The Hamiltonian of the BCi Inozemtsev model is given as follows: 

(2.1) H := -— + J2hik + l)p{x + uJi), 

i=0 

where p{x) is the Weierstrass p-function with periods {2uji,2uj2,), uq = 0, U2 = 
— {uji+ujs) and k {i = 0, 1, 2, 3) are coupling constants. In this paper, we will consider 
eigenvalues and eigenfunctions of the Hamiltonian H. Let f{x) be an eigenfunction 
of H with an eigenvalue E, i.e. 

(2.2) {H - E)f{x) = (^-^ + J2 ^^(^^ + + - /(^) = 0- 

We will make a correspondence between equation ( |2.2| ) and the Heun equation as 
was done in Set z = p{x) and Cj = p{uJi) {i = 1,2,3). By a straightforward 
calculation, we have 

(f , 2 / rf2 1 ^ 1 1 1 \ d 



Hence equation (|2.2|) is equivalent to the equation 



, , \ (P I f I 1 I \ d 1 



ydz"^ 2 \z — ei z — 62 z — e^Jdz A^z — ei){z — e2){z — 62) 



3 



where i' ,i" G {1,2,3}, i' < i" , i' and i" are different from i, /(p(x)) = f{x) and 
C = -E + YlLi + Note that ^ + eg + 63 = 0. 
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Equation ( p.3| ) has four singular points ei, 62, 63 and cxo on the Riemann sphere 
and all singular points are regular. The exponents at z = Ci {i = 1, 2, 3) are (/j + l)/2 
and —k/2, and the exponents at 00 are (/q + l)/2 and 

The Heun equation is a standard form of a Fuchsian differential equation with four 
regular singularities. It is written as 

/ ,N ( f d \ ^ f 1 ^ e \ d aBw — q \ , ^ 

2.4 + - + 7 + 7 ^ + , ,w ^ w = 0, 

\\dw J \w w — I w — tj aw w[w — l)[w — t) I 

with the condition a + /5 + l = 7 + 5 + e. It is easy to transform an arbitrary Fuchsian 
differential equation with four regular singularities into the Heun equation. Thereby 
equation (|2.3| ) is transformed to the Heun equation. 

Conversely, if a Fuchsian differential equation with four regular singularities is 
given, we can transform it into equation (|2.3|) with suitable Cj (i=l,2,3) with the 
condition 61 + 62 + 63 = 0, k {i = 0, 1,2,3), and C by changing the variable z 

and making the transformation / i— > (z — 61)"^ (z — 62)°^ — 63)"^/. Thus the 
relationship between the BCi Inozemtsev model and the Heun equation is made. 

In this paper, we discover properties of the Heun equation by using an expression 
of elliptic functions especially for the case k G Z>o {i = 0,1,2,3). Note that a 
Fuchsian differential equation with four regular singularities with exponents (7i,5j) 
{i = 0, 1,2,3) satisfying 7j — 5j G Z + | is transformed to equation ( p.3| ) with the 
condition k G Z>o. 



3. Bethe Ansatz equation 

3.1. The invariant subspace of doubly periodic functions. In this subsection 
we will find finite-dimensional invariant spaces of doubly periodic functions with 
respect to the Hamiltonian H (see (pTl|)) for the case k G Z>o {i = 0,1,2,3) and 
analyze them. In the terminology of the Heun- type equation (^.31), doubly periodic 
functions in x correspond to algebraic functions in the parameter z{= p{x)). 

Let be the space spanned by meromorphic doubly periodic functions up to signs, 
namely 

(3.1) JF = JF 

ei,e3=±l 

(3.2) JF^,,3 = {fix): meromorphic |/(x + 2uJi) = eif{x), f{x + 2u;3) = 63/(0;)}, 

where {2uji,2u3) are basic periods of elliptic functions. 

Let ki {i = 0,1,2,3) be the rearrangement of k {i = 0,1,2,3) such that /cq > 
ki > k2 > k^{> 0). The dimension of the maximum finite-dimensional subspace in 
JF which is invariant under the action of the Hamiltonian is given as follows. 



HEUN EQUATION 5 

Theorem 3.1. (i) If the number /cq + + ^2 + ^3(= + ^1 + ^2 + ^3) is even and 
^0 + ^3 ^ ki + k2, then the dimension of the maximal finite- dimensional invariant 



subspace in T with respect to the action of the Hamiltonian H (see ( \2. 1\ ) ) is 2kQ + 1. 

(a) If the number ko + ki + k2 + ks is even and ko + k^ < ki + k2, then the dimension 
is ko + ki + k2 - ks + 1. 

(Hi) If the number k^ + ki + k2 + k^ is odd and ko > ki + k2 + k^ + 1, then the 
dimension is 2ko + 1. 

(iv) If the number k^ + ki + k2 + k^ is odd and k^ < ki + k2 + k^ + 1, then the 
dimension zs /cq + + ^2 + ^3 + 2. 

Proof. Firstly, by assuming the existence of the finite-dimensional invariant space W, 
we derive the properties of the space, and later we show its existence. 

Let W he a finite-dimensional invariant space in and T« (i = 0,1,2,3) be 
operators on JF defined by T^*^/(x) = /(— x + 2iUi). From the definitions of the 
operator H and the periodicity of the function which belongs to JF, it follows that 
(T(^))2 = 1, T^T^-'') = T(J')T(^), T«/7 = i7T« and T^^)T^^)t'^^)T^^) = 1 on Set 

(3.3) V= ^ (7^(3))a3(2^(2)^a2(^2.(l)^ai(^j.(0)^)ao . 

ao,ai,a2,a:j£{0,l} 

then the space V is finite-dimensional and invariant under the actions of H and T^*^ 
{i = 0, 1, 2, 3). From the commutativity of the operators T*^*) {i = 0, 1, 2, 3), we have 
the decomposition 



(3.4) V= K 



e0i<:i,<:2,e3 ' 



'!0''!1.E2.E3=±1 

<!0I!1<!2<!3 = 1 



where V,,,,,,,,,,, = {/(x) G V|r«/(x) = ej{x) for z = 0,1,2,3} (eo,ei,e2,e3 = ±1). 

Let /(x) be an element of VeQ^ei,t2,<i3- From the definition of T^*^ and Ko,ei,e2,e3) the 
function f{x + Ui) {i = 0, 1, 2, 3) is odd or even. Hence 

00 

(3.5) fix)= 

n=Mi ^ 

for some Mj G Z. The exponents of the operator H at x = Ui are —k and k + 1. Let 
ai be one of the exponents at x = Ui such that aj + is divisible by 2. 

If c^^^ ^0 for some n such that 2n + < a, then the functions i/-' f(x — uji) 
{j = 0,1,2,...) are linearly independent and there is a contradiction with finite- 
dimensionality. Therefore if 2n + < then c^^^ = 0. 

2nH 

Conversely assume c^*^_^i_,. = for all i(G {0, 1, 2, 3}) and n such that 2n + < 
ftj. Then the function Hf{x) is expressed as Hf{x) = Xlnez^*'*^ ^-h ~ '^i)^"'' 
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0, 1, 2, 3) and c 



2n+- 



for n such that 2n + < ai, because ctj is an 



exponent of the operator H aX x = uji. 

Therefore the space Ko,£i,e2,e3 consists of functions whose coefficients of the expan- 
sion (see ( p.5|) ) at X = ciJj (i G {0,1,2,3}) satisfy ^i-^-- = for all i and n such 



2 < ai. From decomposition ( |3.4| ) and the condition eoeie2e3 = 1, the 



that 2n + 

number + + ^2 + as must be even. 

Now we show that the non-zero finite dimensional vector space Ko,ei,e2,e3 
tto + «i + «2 + ^3 G 2Z<o. Let d 



exists iff 



(3.6) K 



eo,<:i,e2,e3 




aiWN'^l (a,i.)Y^ /<X3{x)^°« 



ct{x) 



p(x)", G Z>o; 

otherwise, 



where ai{x) are the co-sigma functions which are defined in the appendix. Then 
the space Ko,€i,e2,e3 is maximum, i.e. K 



<:0i<:iie2,e3 



C K 



e0)ei:<:2,<:3 



and it is seen that H 



^0,<:i,e2,e3 ^0,ei,e2,e3 • 



Note that the dimension of the space Ko,ei,e2,e3 



is c? + 1. Set 



(3.7) 



V 



«0,<;i,e2,e3 ' 



e0.'^l.'!2.'!3=±l 
'!0<il«2«3 = l 



then the space V is the maximum finite-dimensional space spanned by functions in 
and is invariant under the action of the Hamiltonian H. In particular, we have 

W G V G V . The dimension of the space V is the summation of the dimension of 

the spaces K(,^ei,£2,e3 eoeie2e3 = 1- 

If Iq + h + I2 + h is even, then the dimension of the space spanned by all finite 

dimensional invariant spaces in JF is 



(3.8) 



^0 + ^1 + ^2 + h 



+ 



+ 1 



^0 — ^1 ^ ^2 + ^3 



Iq + h — h — h 



+ 



~ h + h ~ h 



If lo + h + I2 + h is odd, then the dimension of the space spanned by all finite 
dimensional invariant spaces in is 



(3.9) 



lo + h-k + k + l 



-k + li + l2 + h + l 



Hence the theorem is obtained. 



□ 



□ 



Remark Theorem J!T and its proof is a generalization of |T^, §23.41]. 

We will consider the multiplicities of eigenvalues of the Hamiltonian H on the space 
spanned by finite dimensional invariant spaces in JF. We will obtain the following 
theorem, which will be used later. 
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Theorem 3.2. Let V be the space which has appeared in Theorem \3. 1\ (i.e. ( \3. Tj) ). 

If two of li {i = 0, 1,2,3) are zero, then the roots of the characteristic polynomial of 
the operator H (see (\2. 1\) ) on the space V are distinct for generic oji and oj^,. 



Proof. By shifting the variable a; — »• x + cUj and rearranging the periods, it is sufficient 
to show the case I2 = h = Q and /q > k- Set $(2) = {z - eiY^{z - e2T^{z - e3)"^ 
where Ui = -li/2 or {k + l)/2 for each i{e {1,2,3}), and f{z) = <l>{z)f{z). If the 
function f{z) satisfies equation (|2.3|), then 



Note that the exponents of equation (|3.1CI| ) at 2; = 62 are 



(3 10) ^'/(^) ^Y '^^^ + ldfiz) ^ / (ai + a2 + a3-|)(Qi + «2 + a3 + ^)2: 
dz"^ "ill ^ ~ dz I {z - ei){z - e2){z - 63) 

f - ei(a2 + a^y - 62(01 + a^f - 63(01 + 02^] ,, , 

+ - 7 w w ^ f(z) = 0. 

(Z - 6i)(z - 62)(Z - 63) J 

(0, 1/2) (&2 = 0) 

(-1/2,0) (a2 = l/2) • 
Write f(z) = J2'^o (^r(z — 62Y, then the recursive relations for coefficients Cr are given 
as follows: 

(3.11) (62 - ei)(e2 - e3)(2a2 + ^)ci + qco = 0, 

(3.12) (r - 1 + 7i)(r - 1 + 72)c,_i + (62 - ei)(e2 - 63)(r + l)(r + 2^2 + ^)c,+i 
+ (((2a2 + 263 + r)(62 - 61) + (2a2 + 2ai + r)(e2 - 63))r + q)cr = 0, 

where r > 1, 71 = ai + 02 + as - 72 = "i + "2 + "3 + ^ and q = j - 61(02 + 
c^s)^ — 62(01 + 03)^ — 63(01 + 02^ + 627172- Either 71 or 72 is an integer. Let this 
integer be denoted by 7. 

The coefficients are polynomial in E of degree r. We denote by Cr(E). If E is 
a solution to the equation Cr(E) = for the positive integer r = 1 — 7, then it is seen 
that C2-^(E) = from recursive relation (|3.12|) , and also Cr(E) = for r > 2 — 7. 
For the value E, the power series f(z) = Ylr '^r(E)(z — 62Y is a finite summation and 
the function f(p(x)) is meromorphic and doubly periodic. 

If Iq + li is even, it occurs that (ai, 0:2, 03) = (— /i/2, 0, 0), (— /i/2, 1/2, 1/2), ((li + 
l)/2, 1/2,0), ((li + l)/2,0, 1/2) and the degrees of the polynomials Ci^;y(E) are (Iq + 
/i + l)/2, (Iq + Ii)/2, (Iq — Ii)/2, (Zq — ^i)/2 respectively. We denote these polynomials 
by c«(E), c(2)(E), c(3)(E), and c^^^(E). 

Since dj (i = 0, 1,2,3) belongs to |Z, the function $(p(x))/(p(x)) is doubly peri- 
odic up to signs. Let Ei (resp. E2) be a solution to c^^'^\E) = (resp. c^'^'^\E) = 0) 
for some ii (resp. ^2) and let fi(z) (resp. f2(z)) be the corresponding solution to 
equation ( ^.lUI ) which is determined by recursive relations ( |3.11| , |3.12| ). The periodic- 
ity of the functions ^(p(x))fi(p(x)) and ^(p(x))f2(p(x)) , (i.e. the signs of $(p(x + 
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^i))fi{p{x + uji))/<l>{p{x))fi{p{x)) and <l>{p{x + uj,))f2{p{x + u,))/<l>{p{x))f2{p{x)) 
{i = 1, 3)) are different if ii ^ 12- 

Assume that ii^ is a common solution to &^{E) = and c^^\E) = for z 7^ j. 
Let f^^\z) (resp. f^^\z)) be tlie solution to differential equation ( p.lO| ) related to 



the value E. Then the functions /'•*■' (p(a;)) and f^^\p{x)) are the basis of solutions 
to (|2.2|). However, this contradicts the periodicity of the solutions when considering 



exponents. Hence, we proved that the equations c^^\E) = and c^^\E) = (z 7^ j) 
do not have common solutions. 

The degree of the polynomial c^^\E)c^'^\E)c^^\E)c^'^\E) is equal to the dimension 
of the space V (see ( p.7|) ). If we show that the roots of the equation c^''\E) = 
{i = 1, 2, 3, 4) are distinct for generic ooi and uj^, the theorem is proved. 

It is enough to show that if 61,62, 63 are real and {62 — ei)(e2 — 63) < then the 
roots of the equation c«(E) = {i = 1,2,3,4) are real and distinct for the case 
(ai, 62, as) = {-h/2, 0, 0), {-h/2, 1/2, 1/2), {{h + l)/2, 1/2, 0), {{h + l)/2, 0, 1/2). 

We will show that the polynomial Cr{E) (1 < r < 1 — 7) has r real distinct roots 
sf^ {i = l,...,r) such that s^{^ < sf"') < s^^ < s^'^ <■■■< s^jl, < sj^i'^ < s^^^ 
by induction on r. The case r = 1 is trivial. Let k G Z>i and assume that the 
statement is true for r < k. From the assumption of the induction, sf'^ < sf^ < 



s^'' < s''2 < ■ ■ ■ < s^k-i^ < s^k^- It is immediate from ( |3.12| ) that the leading 
term of the polynomial Cr{E) in E is positive. Since sf'^^^ and s^'^^^^ satisfy the 
equation Ck-iisf^'^^) = Ck-i{s\'^^^) = and sf''^^ < sf''' < s\'^^\ the sign of the 
value Cfc_i(4^'') is (— l)'^^*. By recursion relation ( p.l2| ) and the equation Ck{sf''') = 0, 
the sign of the value Ck+iisf'^) is opposite to that of Ck-iisf'^). Combined with the 
asymptotics of Ck+i{E) for E +00 and E —00 and deg^ Ck+i{E) = + 1, it 

follows that the polynomial Ck+i{E) has k + 1 real distinct roots and the inequality 

(fc+i) ^ (fc) ^ (fe+i) / (fc) / / (fc) / (fc+i) . ^- J 
s\ < Si < S2 < S2 < ■ ■ ■ < si ' < sl^-^ IS satisned. 

For the case that /q + h is odd, the proof is similar. □ 

The following proposition was essentially proved in the proof of the previous theo- 
rem. 

Proposition 3.3. Suppose h = h = ^ o-nd lo,h G Z>o- If 1^2 -,63 are real and 
{62 — ei)(e3 — ei) > 0, then the roots of the characteristic polynomial of operator H 
on the space V are real and distinct. 



Proof. It is sufficient to prove for the case Iq > li. If 61,62,63 are real then the 
condition (62 — 61) (63 — 61) > is equivalent to ((62 — 61) (62 — 63) < or (63 — 61) (63 — 
62) < 0). For the case (62 — 6i)(62 — 63) < 0, the proposition is proved in the proof 
of the previous theorem. For the case (63 — 61) (63 — 62) < 0, the proof is similar. □ 
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3.2. Results from the method of monodromy. In this subsection, we present an 
integral expression of solutions to equation ( p.2| ) and related results for the case G 
^>o {i = 0, 1, 2, 3). First, we show a proposition which is related to the monodromy 
of solutions to equation ( p.2| ) for the case k G Z>o {i = 0, 1, 2, 3). 



Proposition 3.4. If lo,li,l2,h £ ^>o then the monodromy matrix of equation ( \2.3( ) 
around each singular point is a unit matrix. 



Proof. Each solution to equation ( p^.2| ) has singular points only at x = nitui + n^tos 
{ni,n^ G Z). Due to periodicity, it is sufficient to consider the case a = uJi {i = 
0,1,2,3). The exponents at the singular point a = Ui are — and /j + 1, hence 
there exist solutions of the form fa,i{x) = (x — a)~''(l + YlT=i'^kx'^'^) and fa,2{x) = 
{x — ay^~^^{l + J2'k'=i'^'k^^'')' because equation ( p. 21) lacks a ffist-order differential 
term and because the functions p{x + uj) (j = 0, 1, 2, 3) are even and because the 
numbers Ij + 1 — {—Ij) are odd. Then it is obvious that the monodromy matrix for 
the functions fa,i{x) and fa,2{x) around the point a; = a is a unit matrix, hence the 



local monodromy matrix of equation (|2.2|) around each singular point is also a unit 



matrix. □ 



We explain how to choose particular solutions to equation (|2.2| ) by applying Propo- 



sition 3.4 



Let Mi (i = 1, 3) be the transformations of solutions to differential equation ( p.2|) 
obtained by the analytic continuation x ^ x + 2uji. Since all the local monodromy 
matrices are unit, the transformations Mj do not depend on the choice of paths. 
From the fact that the fundamental group of the torus is commutative, we have 
M1M3 = M3M1. The operators Mi act on the space of solutions to equation ( p.2|) 
for each E, which is two dimensional. By the commutativity M1M3 = M3M1, there 
exists a joint eigenvector A(a;, E) for the operators Mi and M3. From Proposition 
P^ , the function A(x, E) is single-valued and satisfies equations {H — E)A{x, E) = 0, 
MiA{x,E) = miA{x,E) and M3A{x,E) = m^A{x,E) for some mi, 7713 G C\{0}. We 
immediately obtain relations {H — E)A{—x, E) = 0, MiA{—x, E) = mi^A{—x, E) 
and M3A(-x, E) = m^^A{-x, E). 

Now consider the function S*(a;, E) = A{x, E)A{—x, E). This function is even (i.e. 
E*{x,E) = E*{-x,E)), doubly periodic (i.e. E*{x + 2ui,E) = E*{x + 2uj3,E) = 
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E*{x,E)), and satisfies the equation 



-2 ^ hik + + ^iU j 2*(i, E) = 0, 



which products of any pair of solutions to equation (|2.2| ) satisfy. 

The following proposition is used to determine an appropriate constant multiplica- 
tion of the function 'E*{x, E). 



Proposition 3.5. If lo,li,l2,h G Z>o then the equation 

(3.13) ^^(^^ + + - dx 

-2 (^k{k^l)p'{x^uj,)j^ H(a;,E) = 0, 
has a nonzero doubly periodic solution which has the expansion 

(3.14) S(x, E) = co{E) + fi^M^ + ^^)'"'' 

where the coefficients Cq{E) and b^j\E) are polynomials in E, they do not have com- 
mon divisors, and the polynomial co{E) is monic. Moreover the function S(a;,-E) is 
determined uniquely. Set g = deg^co{E), then deg^b^j\E) < g for all i and j. 

Proof. Recall that the function 'E*{x,E) = A{x, E)A{—x, E) is even and doubly pe- 
riodic, and it satisfies equation ( p.l3|) . 
Here we show the following lemma. 



Lemma 3.6. Suppose /o;^i;^2;^3 £ Z>o and fix E E C If equation l\2.2^ does not 
have solutions in T , then the dimension of the space of even doubly periodic functions 
which satisfy equation ( \3.13i ) is one. 



Proof of the lemma. Since the dimension of even (or odd) functions (not necessarily 
doubly periodic) which satisfy equation ( p2.2| ) is one and solutions to equation ( |3.13| ) 



are spanned by products of two solutions to ( p^.2| ), the dimension of even functions 
which satisfy equation ( |3.13| ) is two. 



Suppose that the function A{x,E) is proportional to A{—x,E), then by definition 
one has mi = m^^ and = m^^, that is ml = mg = 1, and it contradicts to the 
assumption of the lemma. Therefore the functions A(x, E) and A(— x, E) are linearly 
independent. 
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If the dimension of even doubly periodic functions which satisfy equation ( |3.13| ) 
is no less than two, all even solution to ( p.l3| ) must be doubly periodic. Hence the 
function A(a;, E)"^ + A{—x, E)^ is doubly periodic. Since A(a;, E) ^ JF, it follows that 
ml 7^ 1 or TTig 7^ 1. For example, suppose mf ^ 1. From the periodicity of A{x,E), 
A(x, Ef + A(-x, Ey = A(x + 1, Ef + A(-x - 1, E)^ = mlA{x, Ef + m^^A(-x, Ef 
which contradicts the linear independence of K{x,E) and A{—x,E). Therefore the 
dimension of even doubly periodic functions which satisfy equation ( |3.13|) is less than 
two. Since the function i?) is an even doubly periodic function that satisfies 
( p.l3| ), we obtain the lemma. □ 

We continue the proof of Proposition Since the function i?) is an even 
doubly periodic function that satisfies the differential equation ( p.l3|) and the ex- 
ponents of equation ( p.l3| ) aX x = uji are — 2Zj, 1, 2Zj + 2 , it admits the following 
expansion: 

3 h-l ^ 

(3.15) S*(x, E) = do{E) + 

Substituting (|3.15| ) for ( p. 131) , we obtain relations for the coefficients Cq{E) and b^j\E). 
On the relations, the coefficient V~^\e) is expressed as a polynomial in V~^){E) {j' < 
j) and E, and the coefficient Cq{E) is expressed as a polynomial in V~-){E) {i = 
0,1,2,3; 0< J </i-l) and E. 

Through a multiplicative change of scale given by 

3 / — 1 

(3.16) E{x,E) = =co(i?) + EE^f (^)^(^ + ^^)'"'' 

' i=o i=o 

we can choose the coefficients in such a way that Cq{E) and b^j \E) are polynomials 
in E, they do not have common divisors, and the polynomial Cq{E) is monic. 

If there is a function in that satisfies {H — E)f{x) = 0, we have f{x) G Ko,ei,e2,e3 
(see (^.71)) for some eo, ei, €2, ea G {±1}, and the value E is an eigenvalue of H on the 
space Ko,ei,e2,e3- F^om Theorem ^TT] and Lemma any even doubly periodic function 
that satisfies equation (|3.13|) is determined uniquely up to constant multiplication 
except for finitely many E. Hence the polynomials Co{E) and b^j\E) are determined 
uniquely. 

Let E{x,E) = Ylk=o^a'-k{^)'^^ ^ where g' is the maximum of the degrees of cq{E) 
and bf{E) {i = 0,1,2,3; < j < k - 1) in E. Substituting this into (|3l3|) 
and computing the coefficient of i?^'"*"^ yields ^ao(x) = 0, and therefore we have 
'^^gEbf{E) < deg£;Co(^), g' = g, and ao{x) = 1. □ 

We will get an integral representation of A(x, E) by using the function S(a;, E). The 
following argument is motivated by the book of Whittaker and Watson [|T2i §23.7]. 
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Since the functions A(a;, E) and A(— x, E) satisfy differential equation (p.2|), we have 



i,A{x,E)-A{x,E)£ 

p\ _ TP\dK(_^ _ nn(T?\ ^r^A dlogA{x,E) 



£ (A(-x, E)£A{x, E) - A(x, E)±A{-x, E)) = A(-x, E)£,k{x, E)-A{x, E)£a{-x 



0. Therefore A{-x,E)£A{x,E) - A{x,E)j^A{-x,E) = 2C{E) and 
'^^°^^dx^'^^ = 1^^^ for some constant C{E). From the equation ^jx^'^^ di'^'^^ 

1 dS{x,E) , 

^-^^ we have 



z.{x,E) dx 



^3^7^ cilogA(±x,E) _ 1 dEjx^E) ^ C{E) 



(3.18) A(±x,E) = v/2(x,E)exp ^ ^ 



2S(x,£^) dx H(x,i?) 

which is integrated to yield 

± C{E)d x 
(^ 

This formula is the integral representation of the solution to equation 

It remains to determine the constant C{E). Differentiating equation ( |3.17|) yields 

^3^g^ 1 d'Aix,E) ( 1 dA{x,E)^' 



A{x,E) dx^ \A{x,E) dx 

1 d^E{x,E) 1 f 1 d'E{x,E)V C{E) dE{x,E) 



2E{x,E) dx^ 2 \E{x,E) dx J E{x,EY dx 
Set Q{E) = -C{Ef. By substituting (|3) and ([3T7| ) into (|3l9|) , we obtain the 
formula 

(3.20) Q{E) = -C{Ef =S(x, Ef ( £^ - J] k{k + l)p{x + 

\ i=0 

1„, ^,^2^(3;^^) l/dH(x,E)^^ 

+ -^(x,E)- 



2 ' dx^ 4 y dx 

It follows that Q{E) is a monic polynomial in ii^ of degree 2g + 1 from the expression 
for E{x^E) given by ( |3.14| ) and Proposition |3.5| . In summary, we have the following 
proposition. 

Proposition 3.7. Let E{x,E) be the doubly periodic function defined in Proposition 
3.5 and Q{E) be the monic polynomial defined in \3.2(\) . Then the function 

±^-Q{E)dx 
E{x,E) • 



(3.21) A(±x, E) = V5(x, E) exp 



is a solution to the differential equation (\2.^ ). 

From the previous proposition, if is a solution to the equation Q{E) = 0, then 
the corresponding function A(x, E) is an element of JF. 

Conversely, we suppose A(x, E) e J-'. For this case, the function A(— x, E) has 
the same periodicity as A(x,-E'). If the function A{x,E) is neither even nor odd, 
then the function A(x, E) + A(— x, E) is non-zero, even, and doubly periodic up 
to signs, and the non-zero function A{x,E) — A(— x, E) has the same periodicity. 
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The function A{x,E) + A{—x, E) belongs to the space Ko,ei,£2,e3 (see ( p.7|) ) for some 
60,61,62,63 G {±1}- Since the function A{x,E) — A{—x,E) has the same periodic- 
ity and the opposite parity as A{x,E) + A{—x,E), we have A{x,E) — A{—x,E) G 
V^-eg _ei,_e2,-e3- Let ai G + 1} be the number such that aj + is divisible 

by 2, then we have dim V^y^ei,e2,e3 — niax(0, 1 — "o+Qi+"2+"3 ^ g^^^ dim V^.^q -ea = 
max(0, "o+"i+Q2+a3 _ ]^^_ Hence either the space Ko,ei,£2,e3 the space Vleg^-e^ 
is {0}, which contradicts the property A{x, E) 7^ ±A(— x, -E). Therefore the function 
A{x,E) must be odd or even. If the function A{x,E) is odd or even, the functions 
A(a;, E) and A(— x, E) are linearly dependent. Thus the Wronskian A(— x, E)-^A{x, E)- 
A(x, E)£A(-x, E) is equal to and C{E) = Q{E) = 0. 
Therefore the following theorem is proved: 



Theorem 3.8. The equation has a solution in the space T (^^) if and only if 
the value E satisfies the equation Q{E) = 0. 

Let us count the degree of the polynomial Q{E) and discuss the relationship to the 
maximum finite-dimensional invariant subspace V (see ( |3.71 )) in JF. 



Conjecture 1. (i) The degree of the polynomial Q{E) in E is equal to the dimension 
of the space V. (See ( ^^) and Theorem \3. J\ ) 

(a) For generic Ui and u^, roots of the equation Q{E) = are distinct. 

Proposition 3.9. Conjecture [J is true if two of li {i = 0, 1, 2, 3) are zero. 

Proof. It is sufficient to show the case ^2 = ^3 = and Iq > h. 



The function S(x, E) admits expression (|3.14|), but we now use the other expression 



-| lo+h 

(3.22) E{x,E) = _ Yl ^jm^-eiY, 

where z = p{x). We are going to obtain a recursive relation for the coefficients aj{E). 

Since E{x,E) satisfies equation (|3.13|) , we find that 
(3.23) 

j(j - 2/0 -l)(2j- 2/0 -l)a,(^) 

+ 2(j - /o - l)i~E + (Sf - 6(/o + l)j + 2/2 + 5/0 + 3)ei)a,„i(E) 

+ (2j - 2/0 - 3) (j - 2 - /o - - 1 - /o + h){ei - 62) (ei - e3)a,_2(^) = 0, 

where j = 1, 2, . . . , Zq + + 1 and a_i(_E) = aig^i^^i{E) = 0. 

By recursive relation ( ^.23| ), we have 
(a) All coefficients aj{E) {j = 1, . . . Zq + ^i) are divisible by ao{E) as a polynomial in 
E, 61,62, 63. Hence we may set aQ^E) = 1. 
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(b) The coefficient aj{E) (j = 0, 1, . . . , lo + h) is a polynomial in variables E, ei, 62, 63 
and the homogeneous degree of aj{E) w.r.t. E, ei, 62, 63 is j. 

(c) deg^aj(E) < j. 

(d) If J < lo then deg^ajiE) = j. 

Put j = lo + h + l in ( ^.23| ), then deg£;a,„+i,_i(^) < deg^ ai,,+i^{E) + 1. By putting 
j = lo + li — k + 1 (fc = 1, . . . , Zi — 1) in ( |3.23| ), we also obtain deg^ ai^+i-^^k-iiE) < 
degg a;Q+/j_fc(-E) + 1. It follows that deg^ (ij{E) < Iq for all j, deg^ 0'j{E) = J = 
Iq, and the coefficient of E'-° in ai^^{E) does not depend on ei, 62, 63. 

From relation ( p.20| ), it is seen that deg^ Q{E) = 2/o + 1 and the coefficient of 
^2io+i jj^ Q{E) does not depend on ei, 62, 63. From Theorem |3.8| , equation (|2.2|) has 



a solution in JF iff the value E satisfies the equation Q{E) = 0, and from Theorem 
|3.1| , dimV^ = 2/0 + 1 (see ( p.7|) ). From Theorem p.2| , the roots of the characteristic 
polynomial of the operator H (see ( p.l|) ) on the space V are distinct for generic ui and 
uj^. Hence the characteristic polynomial of the operator H on the space V coincides 
with the polynomial Q{E) up to constant multiplication. 

Therefore we obtain the proposition. □ 

The following proposition is proved case by case. 

Proposition 3.10. Conjecture [| is true for the case lo + h + h + h ^ 8- 

3.3. Bethe Ansatz equation in terms of sigma functions. In the previous sub- 
section, we introduced the elliptic function S(a;,i?) and the polynomial Q{E). In 
this subsection, we will propose the Bethe Ansatz method. We postulate the Ansatz 
(assumption) that the eigenfunction of the Hamiltonian ( p.2| ) is written as ( |3.31| ) 
with auxiliary variables [ti, . . . ,ti). The Bethe Ansatz method replaces the spectral 
problem (see ( |2.2| )) with solving transcendental equations ( p.29| ) for a finite number 
of variables (ti, . . . ,t;), which we call the Bethe Ansatz equation. It will be shown 
that the basis of solutions to equation (|2.2|) is written as the form satisfying Ansatz 



( p.31| ) except for finitely many eigenvalues E. The precise statement will be given in 



Theorem |3.12| . Note that the case of the Lame equation (i.e. the case h = h = h = 0) 



was treated in [[1^] 



Throughout this subsection, we assume Q{E) 7^ 0. As was discussed in section |3l^ , 
the functions A(a;, E) and A{—x, E) are linearly independent. 

For simplicity, assume Iq 7^ 0. The function S(x, E) is an even doubly periodic 
function which may have poles of deg X = cUj (i = 0, 1, 2, 3). Therefore we 

have the following expression; 



for some values ti, 
depend on E. 
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,ti 



Proposition 3.11. Assume Q{E) 7^ 0, then the values tj, —tj, Ui (j = 1, . . . i 



0, 1, 2, 3, see (|3.24| )) are mutually distinct up to periods of elliptic functions. 



Proof, (a) We show tj 7^ uji up to periods of elliptic functions for each j G {1, . . . , /} 
and i G {0,1,2,3}. 

Since the function E) is an even doubly periodic function and satisfies equation 
( P.13| ), it has a pole of degree 2/i or a zero of degree 2/j + 2 at x = ui {i = 0, 1, 2, 3). 
Suppose that the function S(a;,i?) has a zero of degree 2/j + 2 at x = Ui for some 
i G {0, 1, 2, 3}. Then the functions A(x, E) and A(— x, E) have zeros of degree k + 1 
at X = Ui, because E{x,E) is proportional to the product A{x, E)A{—x, E), A{x,E) 
and A(— X, E) satisfy differential equation ( p.2| ) and their exponents at x = tUj are 
— /j and li + 1. From the condition Q{E) 7^ 0, the functions A{x,E) and A(— x,-E) 
are linearly independent and form a basis of solutions to equation (|2.2| ). It follows 



that every solution to equation (|2.2|) has a zero of degree k + 1 at x = Ui and this 
contradicts the fact that one of the exponents at x = cuj is —li. Therefore the function 
S(x, E) has a pole of degree 2/j at x = Ui {i = 0, 1, 2, 3). 

If = uJi or — ciJj for some j and i, then from expression ( |3.24|) it contradicts the 
degree of the pole at x = t^j. 

(b) Now we show tj ^ tj, {] ^ f, G {1, . . . , /}) and t^ ^ -t^, (j, j' G {1, . . . , /}) 
up to the periods of elliptic functions. If tj = —tj, then we have tj = Ui for some 
i G {0, 1, 2, 3} from the condition 2tj = mod 2u;iZ + and it is reduced to the 
case (a). 

Assume {tj = tj> {j ^ f, G {!,...,/}) or tj = -tj, (j,/ G {1, ...,/})) 

and tj 7^ oji (j G {1, . . . G {0, 1,2,3}), then the function H(x, -E) has zeros of 
degree no less than 2 at two different points x = tj and x = —tj. The degree of 
zeros of the functions A(x, E) and A(— x, E) at points x = tj and x = — are or 1, 
because A(x, E) and A(— x, E) satisfy differential equation (|2.2|) which is holomorphic 
at the points x = ±tj. Since H(x,-E) is proportional to A{x, E)A{—x, E), there 
is a contradiction with the degree of zeros of the function S(x, -E). Therefore the 
proposition is proved. □ 



Set z = p(x) and Zj = p(tj). From formula (|3.20|) , we obtain 



dE{x,E)y ACiEf 



|2 • 



16 



KOUICHI TAKEMURA 



We fix the signs of tj by taking 



dE(x,E)\ 2C{E) 



If we put 2C{E)/'E{x, E) into partial fractions, it is seen that 



2C{E) 



where ({x) is the Weierstrass zeta function. It follows from formula ( p.l8|) that 



(3.25) A{x,E) 



{p{x) - ei)'i(p(x) - e2)'2(p(x) - es)'^ 

-x)-2xC(t,)) 



exp I - ^ (log a{tj +x)- log a{tj 



; exp -X > Qitj. 

a(x)'oai(x)^ia2(x)'^a3(x)'3n' cr(t,) j^,'" 



where a{x) is the Weierstrass sigma function and (Tj(x) (i = 1,2, 3) are the Weierstrass 
co-sigma functions. Consider fixing the value E. If A(x, E) ^ JF, then A(x, -E) 
and A(— a;, E) are linearly independent and the linear combinations of A(x, E) and 
A(— X, i?) exhaust solutions to equation (pl2|). Hence if there is no solution to ( p.2|) in 
JF, a basis of the solution is written as A(x, E) and A(— x, E) where A(x, E) is given 
as ( p.25|) . Note that the condition for E that there is no solution to (|2.2|) in T is 
equivalent to the condition that the value E satisfies Q{E) ^ 0. 

Conversely assume that differential equation (|2.2|) has the solution written as 



^'■''^ = a(x)^o"r,)"l^(^)t3(,)^3 ^-p(^^) 
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for some c and tj (j = 1, . . . , /) such that tj ^ t,y (j 7^ j') and tj 7^ (j = 1, . . . , /, i = 
0, 1, 2, 3). By a straightforward calculation, it follows that 

(3.27) 



A(x) I ^ k{k + + oji) + 2/oC(a;) -c - ((^i) 

. i=0 \ j=l 



+ 2 ^ liC{x + tUi) ( -c - /C('^i) - X] C(^i - 
i=i \ j=i 

+ 2j2cix+ tj) ( c + ^^^^ - ^^■) - ^oC(-t.) - 5^ hici-tj + uj,) - c(^^)) 

j=l \ k^j i=l 

I 3 

J=l i=0 j<k 

- {loh + hhYi - {loh + ^1^3)62 - (/o^3 + ^1^2)63 + ^ hViC^c + 

i=l 

Note that we used the formula 

(3.28) {({x + a) - C{x + b)- ({a - b)^ = p{x + a) + p{x + b) + p{a - b) 

to derive equation ( |3.27| ). 

Hence we find that the function A(a;) satisfies equation ( ^.2[ ) if and only if tj {j - 
1, . . . , /) and c satisfy the following equations; 

3 

(3.29) J2 + ^'^) - ^oC(-^i) - Yl + - ^(^')) = (j = 1, . . . , /), 
(l-Vo)(c + $^C(t,))=0, 

(1 - 5i^,o) (c + /C(^^) + Yl + ^^■)) = 0' = 1' 2' 3)- 

The eigenvalue E is given by 

3 

(3.30) E=-c^ + {loh + l2h)ei + {loh + hh)e2 + (/0/3 + ^1^2)63 - lM2c + lr]i) 

i=l 

I 3 

j=l i=0 j<k 

Note that there are other expressions of E in terms of {j = 1, . . . ,1) and c. 
Summarizing, we obtain the following theorem: 
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Theorem 3.12. Set 

117=1 + 

(3.31) X{x) = —— — ^--^ —- --^exp(ca;), {l = lo + li + l2 + h) 

and assume tj ^ tji {j ^ j') and tj ^ uJi {j = 1, . . . ,1, i = 0,1, 2, 3). 

The condition that the function A{x) satisfies equation ( ^.Sj ) for some E is equiva- 
lent to tj {j = 1, . . . ,1) and c satisfying relations 2i\ ), with the eigenvalue E given 
as ( 'J. 3(\ ). For the eigenvalue E such that Q{E) 0, a basis of solutions to l ^.^ 



is written as A(x) and A{—x), and tj (j = 1,. . . ,1) such that tj ^ tj/ (j 7^ j') and 
tj T^uj, {j = l,...,l, i = 0, 1,2,3). 

Remark (i) Comparing with 0, 3, it would be reasonable to call the function A(x) 
the Bethe vector and equation ( |3.29|) the Bethe Ansatz equation. 



(ii) For each Iq, li, I2, 13, the number of eigenvalues E which do not have a basis of 
solutions to (|2.2|) as the Bethe vectors (p.31|) is at most the dimension of the space V 



(see ( p.Tp ) which is given in Theorem |3.1| . 

(iii) If 7^ and h then the number of equations ( |3.29| ) is strictly greater than 
the number of the variables tj (j = 1, . . . , /), although equations ( |3.29| ) have solutions 

{t,,...,ti) ifg(E)^o. 

For the case h = h = h = ^ (the case of the Lame equation), the equation 
c + X]j=i C{tj) = is obtained by summing up other equations in ( p. 291 ). Hence the 
number of equations is essentially equal to that of the variables tj {j = 1, . . . ,Iq). □ 

3.4. Bethe Ansatz equation in terms of theta functions. We will introduce 
the Bethe Ansatz equation written in terms of theta functions. 

Set T = UJ3/UJ1 and Ui = 1/2, then the Hamiltonian (see ( p.lj) ) is described as 

(3.32) H = -^ + lo{lo + l)p{x) 

+ hih + l)p + + ^2(^2 + l)p (^x + + /3(/3 + l)p [x + ^ 

By replacing the sigma function with the theta function by formula ( [A.7| ), it is seen 
that if Q{E) 7^ 0, then the basis of solutions to ( |2.2| ) is written as A{x) and A{—x), 
where 

exp(7r-v/^cx) n',- T 0(x + tj) 

^ ' ^ ' e{xy'^e{x + i/2y^d{x + {i + T)/2y^e{x + T/2y^' 

tj 7^ tj, {j ^ f), tj ^ iZ + f Z (j = 1, . . . , and l = lo + h + l2 + h. 

By similar calculations as the previous subsection, we obtain relations for tj {j = 
and c for which the function A{x) becomes an eigenfunction of operator 
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Theorem 3.13. The function A(x) (see ^3. 3^ ) with the condition tj ^ ty (j ^ j') 
and tj ^ |Z + is an eigenfunction of ^3.32[ ) if and only if tj (j = 1, . . . , l{= 
Iq + h + h + h)) (I'lT'd c satisfy the relations, 



(3.34) 



CTT^ 



- , , , e'jt,) , , e'jt, - 1/2) 



eit,-i/2) 

e'{t,-{l + T)/2) e'{t,-r/2) 

+ '3 



0, (j 



^(t,-(l + r)/2) ■'eit,-T/2) 
(1 - Ko) (^cnV^ + vr v^(/, + ^s) + ^ ^] = 0' 

V V 12 3; Z.^(t^._i/2) 



-0 



(1 - 5h,o) 



(1 - 5j,,o) CTTV^ - 7rv^(/o + ^1) + 5^ 



g-(t,-(l + r)/2) 
. e(t,-(\^T)l2) 



0. 



The eigenvalue E is given by 



(3.35) E = 7i'{c^ + (/o + /i)(/2 + I3)) + /(/ + l)Vi - Yl 



j<k 



0"(tj — tk) 
9{tj — tk) 



+ (Wi + ^2^3)ei + {loh + kh)e2 + {kh + /i/2)e3 



.7=1 



"{t,) , , 9"{t,-l/2) , , d"{t,-{l + T)/2) , , e"(t,-r/2) 



+ /3 



e{t,) 'e{t,-ii2) ' e{t,-{i + r)/2) "eit,-T/2) 



If the value E satisfies Q{E) 7^ 0, then a basis of solutions to ^2.3i) is written as 
A(x) and A{—x) , where the function A{x) is given in 33i ) and satisfies the condition 
tjT^tj' (jV/) andt^^^Z + ^Z (j = l,. ..,/). 



It would be hopeless to attempt to solve equations ( |3.34D explicitly for generic r. 
Instead, it may be possible to solve them for the case r ^/—loo and to gain some 
knowledge for the case r 7^ V— loo from this limit. We will investigate solutions to 
equations (p.34|) for the trigonometric (r — > \^—loo) case in the next section. 



4. Trigonometric limit 



4.1. Trigonometric limit of the polynomial Q{E). In this section, we will con- 
sider the trigonometric limit r V— loo. 
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The trigonometric limits of the elhptic functions are given by 
(4.1) p{x)^-— + — -, p{x+- 



3 (sinvrx)^ \ 2j 3 (cosTrx)^' 

T\ 71^ [ 1 + t\ tt^ e{x) 



(4.2) + ^-sin^rx, 

as r ^ V— loo. They converge uniformly on compact sets. Due to this limit, the 
term ^2(^2 + l)p(a^ + ^^) + ^3(^3 + + f ) converges to a constant. From now 

on, we consider the case h = h = 0. Set p = exp(27rA/— Ir). Note that the limit 
r V— loo corresponds to the limit p 0. 
As p — > 0, we obtain H Ht — Ct, where 

j2 2 2 

(4.3) Ht = + /o(/o + 1) ,."" ,2 + ^1(^1 + 1) / X2 > 

dx"^ (smvrx)"^ (cosvrxj^ 

(4.4) C'T = (/o(/o + l) + /i(/i + l))7rV3. 

We will consider the trigonometric limit of functions which have appeared in section 
|3^ . Note that ei 62 — > and 63 — > as p ^ 0. The function S(x,-E) 

(see ( |3.14|) ) admits the trigonometric limit 

(4.5) S,(..£)=c„(B) + 5:— ^-y-y + E 



,sm7rx)^('o ■'^ (cos7rx)^('i 
i=o ^ ^ i=o ^ ^ 



and the function Sr(x, i?) satisfies the trigonometric version of equations ( |3.13| ) and 
( ^.2U| ). We will calculate the polynomial Qt{E) which is the trigonometric version 
of Q{E). From the trigonometric version of equation ( |3.2CI| ), we have Qt{E) = 
{E + CrycoiEY. 

Proposition 4.1. If Iq + li is even, then 

ip+h I'o-'il 

2 2 

(4.6) Qt{E) = b-iE + CT)ll{E- {2if + Crf \{ {E - {2i - if + Ct)\ 



i=l i=l 



where Ct is defined in jj-A ) CL^d h is a constant which is independent of E. 
If Iq + li is odd, then 

I'o-'il-i 'o+'i+i 

(4.7) QT{E)=b-{E + CT) n {E-{2tf + CTf H (E - {2t - if + Ct)' . 

i=l i=l 

Proof. For simplicity, we prove for the case /q > h and /q, ^1 ^ 2Z>o. The proofs for 
the other cases are similar. 

From the proof of Proposition ^ , it follows that the polynomial Qt{E) is obtained 
by a suitable limit of the characteristic polynomial of the space V (see ( ^.7| )). For 
this case, we have V = Vi, 1,1,1 © VI, 1,-1,-1 © VI, -1,-1,1 © V^,-i,i,-i- 
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Set e = (e2 — 63), then e ^ as r ^ V— loo. We will solve recursive relation 
( p.l2| ) for the case e — > 0. If e 7^ then the invariant subspace Vi, 1,1,1 is given as 

'O+'l f fi-A"'^ 

^1,1,1,1 — ©n=o ^ \ ^^) Then the characteristic polynomial of the space 

^1,1,1,1 ^ ~^ is expressed as 

'o+'i 

(4.8) e-^-Vi n - (20' + Ct) + ^ + . . . , 

i=0 

where C\ is a constant. Similarly, the leading terms of the characteristic polynomials 

of the spaces Vi,i,_i,_i, Vi,_i,_i,i, Vi, -1,1,-1 as e are given as 

(4.9) 

^0+^1 ^Q-^l ^Q-^l 

2 2 2 

n (i? - (20' + , n (i?-(2^-i)'+CT), n (i?-(2^-i)'+CT), 

1=1 i=l j=l 

up to constant multiplications respectively. 

By multiplying the leading terms of these four cases, we obtain the proposition. □ 

4.2. Bethe Ansatz equation for the trigonometric model. In this subsection, 
we will obtain the trigonometric version of equations ( p.34| ). 

Set X = exp(27rV— la^) and Tj = exp{2Tr\^—ltj). Now we state the trigonometric 
version of Theorem |3.1tj| . Set 

and assume Tj 7^ Tj/ (j 7^ j') and Tj ^ 0, ±1 (j G {1, . . . ,lo + h})- Suppose that Tj 
and c satisfy relations, 

(4.11) c + 5^^±^ + /o^ + /i^ = (j = l,...,/o + /i) 

(4.12) (l-5.o,o) =0, {1-S,„o)[c+J2y 



0, 



j=l 1 J \ j=\ 3 



then it is shown that the function Aj'(e ) is an eigenfunction of operator Ht 
([4.3|) with an eigenvalue E = vr^c^. 



Conversely, if the function At(X) = (7riV„(x+iy! (^^- ^ T,, (j ^ j') and 
^ 0, ±1 (j = 1, . . . , /o + ^i)) satisfies the equation i7rAr(e2'^^^) = EAT(e2^^^), 
we obtain equations ( [4.11| , [4.12| ) and = vr^c^. 

Proposition 4.2. Fzx i/ie value E. If Qt{E) 7^ then the basis of the solutions to 
the equation HTf{x) = {E + CT)f{x) is written as Axie'^'^^^^) and Aj'(e~^''^^'^') 
where Aj'{e'^'^^^^) is given in ( U-iQ) with the conditions = '^^fi , Tj 7^ Tji {j 7^ j') 
andT,y^O,±l (j e {l,...,lo + h})- 
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Proof. From expression ( |4.5| ), it follows that 

Co(-E')(sin7rx)^^'"''"'^^ + (lower terms in (sinvrx)^) 



(4.13) Et{x,E) - / • N9; / ' 

(sm7rx)2'o(cos7ra;)2'i 

for some ti, . . . , tjg+j^, where the polynomial Co{E) is common in ( [4.5[ ) and ( [4.13[ ). 
From the relation Qt{E) = {E + Ct)cq{EY , we have Co{E) ^ 0. Hence Sr(x, i?) is 
written as 

5o(^)nS^((smvrx)2-(sinvrt,)^) 



(4.14) Et{x,E) 



(sin TTx)'^'-° (cos 7ra;)^'i 



we 



(4-15) At(x) = exp(7rv^c), 



for some ti, . . . ,i«o+'i- -^y a similar argument to that performed in section 
obtain 

Lfc=i 

(^sm 'n'xy°[cos ttx^ 

where the value c satisfies = -^^r^ and the signs of ti, . . . , t^g+jj are suitably chosen. 

If At{x) satisfies the differential equation {Ht — Ct — E)At{x) = 0, At{—x) also 
satisfies it. From the trigonometric version of formula ( p.l8|) , it is seen that the 
functions At{x) and At{—x) are linearly independent. By a similar argument to the 
proof of Proposition p.ll| , the values ±tj, (j = 1, . . . ,/o + ^i), 0, 1/2 are mutually 
different. By setting X = exp{27i\/^lx) and Tj = exp{27i\^—ltj), we obtain the 
proposition. □ 

4.3. Solutions to the Bethe Ansatz equation. We will solve the trigonometric 
Bethe Ansatz equation. 

Theorem 4.3. Assume Qt{t^'^c^ -Ct) ^ 0. Let {T^, . . . ,Tj^) {I = lo + h) be a solution 
to equations U-^W o-nd let cTj = X]ji<...<j ^j^i • • • -^j^ ^"^^ elementary 

symmetric function of the solution (T°, . . . ,T°). Then there exists a solution to l \4-l^ , 
uniquely up to permutation of (T°, . . . ,T°). The functions ai are determined 



by the following recursive relation, 
(4.16) 

{h - li){c - Iq - h) 
^0 = 1, ^1 = , 

(4.17) 

(/o -li){c + 2i-2-k-h) {i-2~lo-h){c + i-2-lo- h) 
(Ti = .,. . cTj-i H ai-2 

t[t + C) t[t + c) 

/or z = 2,3, ... ,/o + ^1- 

Proof. By the assumption Qt{tt^c^ — Ct) 7^ and Proposition [4.2| , there exists a solu- 

(x-iyo{x+i) 



tion to the equation {HT-n^c^)ATie^^'^'') = of the form At (A) = 
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Then the values (Ti, . . . ,Ti) satisfy equations ( [4.11| , [4.12| ) and we have 
(4.18) 



(X-l)'o(X + l)'i 

By comparing the coefficients of (j = 1, . . . , Z) of the equation {Ht—h'^c'^) At{X) = 
{X = e^'^^^^), we obtain recursive relation ( |4.17 ) for 0",. □ 



Now we exhibit explicit expressions of am- By straightforward calculations, we 
obtain the following proposition: 



Proposition 4.4. (i) If li = then 

m—l 



(4.19) 

(a) If li = 1 then 
(4.20) am = 



a„ 



n 

j=0 



{h-3){c + ]-h) 
(j + l)(c + j + l) 



nr=o'(^o-j)nr=f(c+j-/o) 



(c - /o - l)((/o - 2m + l)c - /q - ^0 + 2m/o + 2m - 2m^) 



(Hi) If li = 2 then 
(4.21) am 



(c-/o)(c-/o-2)- 



(c2(Z2 - (4m - 3)/o + 4m2 - 8m + 2) 

+ 2c{-ll + (4m - 3)/o - (6m2 - 10m + 2)/o + 4m^ - 12m^ + 8m) 

+ 1^- (4m - 3)/j^ + (8m2 - 12m + 1)/^ 

- (8m^ - 20m2 + 8m + 3)/o + 4m^ - 16m^ + IGm^ - 2). 



(iv) If li = Iq then 

(4.22) a2m~i = 0, as, 

(v) If li = Iq — 1 then 

(4.23) a2m-i 



m—l 

n 

j=0 



(j-/o)(c + 2j-2/o) 
(j + l)(c + 2j + 2) ' 



nr=/(j-^o)nr=V(c+2j +1-2/0) 



rm— 1 



(^-i)!nr=o'(c+2j + i) 
nr=iO'-^o)nr=o'(^+2j + i-2/o) 
^!nr=o'(c+2j + i) 
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(vi) If li = Iq — 2 then 

2 UT=2ij - ^o) nr=o'(c + 2j + 2 - 2/o) 

(4.24) a2m-i = 



^-l)!(c+l)nr=o'(c + 2j + 2) 

nr=2(j-^o)nr=o'(^+2j +2-2/0: 

^2m — ^m-l 



^!(c+l)nr=o(c + 2j + 2) 
(c(/o - 2m - 1) + (4m + l)/o - (2m + l)^). 



Remark Proposition |4.4| (i) reproduces Proposition 4.1 in 1^]. □ 

5. ElGENSTATES FOR THE BCi INOZEMCHEV SYSTEM 

5.1. Eigenstates for the trigonometric Hamiltonian. Let us find square-integrable 
eigenstates of the Hamiltonian Ht- Set 

(5.1) $(x) = (sin7rx)'°+^(cos7rx)'^+\ Ht = ^ix)-^HT^x), 

then the gauge transformed Hamiltonian is written as 

d"^ /(/o + 1) cosvrx (/i + 1) sin7rx\ d 22 



(5.2) Ht = --r^ - 27r '-^^—^ '-^-^ + + + 2)^7r^ 

dx"^ \ smvrx cosvrx / dx 

By a straightforward calculation we obtain 

(5.3) Ht(cos 27rx)'" = 7r2(2m + k + + 2)2(cos 2'kx)'^ 

+ 47r2((-m(m - 1) + m(/o - /i))(cos 27rx)""^ 

Therefore the operator TIt acts triangularly on the space spanned by {1, cos 2tix, (cos 27rx)^, . . . }. 
Since (2m + /q + ^1 + 2)^ 7^ (2m' + /q + ^1 + 2)^ (m, m' G Z>o, m 7^ m'), the operator 
TIt is diagonalizable and 

(5.4) HT^mix) = TT\2m + I0 + I1 + 2)2^^(x), 

for some functions iprn{x) = (cos 27rx)™'+ (lower terms). Set w = {1 — cos27ra;)/2. By 

a change of variable, it follows that 

(5.5) 

Ht - 7r2(2m + k + h + 2f = -An^- 

(^2/2/13 \ d 

wil-w)^+ ((/o + /i + 2) + l)w;j — + m(m + /o + /i + 2) 

Hence the operator Ht — vr^(2m + lo + h + 2)^ is transformed to the hypergeometric 
operator. From expression (|5.5| ), we obtain 

, / \ ~ ^ fi , 2Zo + 3 l-cos27rx' 

(5.6) ^rn{x) = CraGm Uo + h + ^, ^ ; ^ 

where the function Gm(«, P', w) =2 Fi(— m, a + m; P; w) is the Jacobi polynomial of 
degree m and is a constant. 
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We define the inner product by 



(5.7) {f,g) = / dxf{x)g{x), 

Jo 

then the values {^{x)ipm{x),^{x)ipm'{x)) are finite and 

(5.8) n\2m + lo + h + 2)2($(a;)7/'„(x), $(x)^„,(x)) 

= {HT^x)^m{x),<^{x)■^m'{x)) = {^x)^|Jm{x) , HT^x)^m' (x)) 

= 7i\2m' + I0 + I1 + 2)2($(x)V^„(x), <l>ix)iJm'{x)). 

Therefore the functions ^{x)ilJrn{x) {m G Z>o) are orthogonal. 

Remark (i) The values and {^{x)iljrn{x),^{x)iprn{x)) are known exphcitly. 

(ii) The set of functions ^{x)ipm{x) {ni G Z>o) is a complete orthogonal basis of the 
space L2(M/Z). For the proof, see § or [0]. 

(iii) The functions ipm{x) are simply the SCi-Jacobi polynomials. □ 

5.2. Jacobi polynomial, Bethe Ansatz and perturbation. In this subsection, 
we will see the relationship between the Jacobi polynomial and the function ob- 
tained by the Bethe Ansatz method. After that, we will show the holomorphy of the 
perturbation with respect to the parameter p = exp(27rV— It) ~ for the Hamil- 
tonian of the Inozemtsev model ( p.l|) with the condition h = h = 0. For the case 
^1 = ^2 = ^3 = 0, the corresponding result is obtained in p. 

Let us recall the Bethe Ansatz method for the trigonometric model, which was 
discussed in section |[ 

Let m be an non-negative integer, and set c = /o + ^i + 2 + 2m, then Qt{t^'^c^ — Ct) 7^ 
0, where Qt{E) and Ct were defined in section Hence the solution Ti, . . . , Ti^^i^ 
of the Bethe Ansatz equation (see ( [4.11| , |4.12| )) with the condition c = /g + ^1 + 2 + 2m 



exists uniquely up to permutation. In this setting, the function 

smn{x + t,) 



rlo+l 

(5.9) At{x^- J^-v^cx 



(sin'7rx)'o(cos'7rx)'i ' 

satisfies HtAt{x) = 7r^c^AT'(x), where exp(27rV— l^j) = Tj and the values (Ti, . . . , Ti^+i-^) 
are a solution to equations ([4.11| , [4.12| ). From the condition QT{7!''^{lo + h + 2 + 2m)^ — 



Ct) 7^ 0, it follows that the functions At{x) and At{—x) are linearly independent. 
Set Af^'^ix) = At{x) - (-I)'oAt(-x), then Afy'^ix) and the degree of the pole 
of the function Afff"^{x) aX x = n (resp. x = n + 1/2) {n G Z) is less than Iq (resp. 
li). Since the exponents of the differential equation {Ht — 7r^c^)A^^™(a;) = at 
X = n (resp. x = n + 1/2) are —Iq and /q + 1 (resp. —li and li + 1), the function 
Aj^^{x) has zeros of degree Iq + 1 (resp. li + 1) sX x = n (resp. x = n + 1/2). 
Since the eigenvalue of the function $(a;)^/'m(a;) with respect to the Hamiltonian Ht 
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is 7r^(/o + /i + 2 + 2m)^ and the function ^{x)ilJrn{x) is holomorphic at x = 0, it follows 
that Afp"^{x) = Bm^{x)ipm{x) for some non-zero constant B^- 

Now, we will show that for each c = Zq + + 2 + 2m {m G Z>o) there exists p G M>o 
such that the solution to Bethe Ansatz equation ( |3.34| ) exists for each p {\p\ < p) and 



it converges analytically to the solution to the Bethe Ansatz equation for the case 
p = 0. Set Cm = lo + li + 2 + 2m. Since the coefficients of the polynomial Q{E) defined 



in section ^]2| are holomorphic in p, there exists ei,pi G ]R>o such that Q{E) ^ for 
all E and p such that \E + Ct — tt'^c^I < ei and \p\ < pi. Let us recall that the 
function E) is expressed as 

_, „ nyi''(p(^)-pfe)) 

Set ttj = p(tj) and z = p{x), then the values aj depend on the values E and p. 

The numerator of equation ( |5.1CI| ) is written as a polynomial in z and the coefficients 
of z^ (j = 1, . . . , /q + /i) are holomorphic in E and p. Since the zeros of the numerator 
of ( |5.10| ) in z are simple, the values aj are holomorphic in E and p such that |i? + 
Cr-vr^c^l < ei, IpI <pi. 

Suppose |£' + Ct — vr^c^l < ei and \p\ < pi, then Q{E) ^ and there exists a solu- 
tion to Bethe Ansatz equation (|3.34|) . Denote the solution by {ti{E,p), . . . , ti^^^i^{E ^p)). 
If we choose the index suitably, then = p{tj{E,p)). 

From Proposition |3.11| and the proof of Proposition [4.2| , it follows that the values 
tj{E,p) are mutually distinct and are not contained in < _|_ (p 7^ 0) 

E and p such that \E + Ct — tt^c^I < ei, \p\ < pi. Since tj(E,p) ^ |Z + |rZ, we 
have p'(tj{E,p)) 7^ and the values tj{E,p) are holomorphic in and p because the 
values aj are holomorphic in E and p. From (|3.35| ), the value c satisfies the relation 
TT^c^ = h{E,p) where 

(5.11) 



h{E,p) = E-{lo + h){lo + h + l)r]i - lohei 

lo+h 



^ - h{E,p)) ( e"{t,{E,p)) e"{t,{E,p) - l/2) \ 

^ 9{t,{E,p)-t,{E,p)) 2^^['9(t,{E,p)) + '9it,{E,p)-l/2)J- 



We are going to investigate the Bethe Ansatz equation with the condition c = Cm 
for each p {\p\: sufficiently small). For the case p = 0, we have h{E — Ct, 0) = E and 
H(7r^c^ — Ct, 0) = 1. By the implicit function theorem, there exists p G M>o such 
that there exists a holomorphic function E{p) in p on the domain \p\ < p such that 
h{E{p),p) = TT^cl and E{0) = tt'cI - Ct. Set i,ip) = tj{E{p),p) (j = 1, . . . , Zq + /i), 
then the function tj{p) is holomorphic in the variable p on \p\ < p, (ti(p), . . . , ^io+h (p)) 
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is a solution to Bethe Ansatz equation ( |3.34| ) with the condition c = c^, and the 
corresponding eigenvalue E{p) is also holomorphic in p. 

Hence we obtain that the solution to Bethe Ansatz equation ( |3.34| ) at p{= exp(27r-\/— 1 
is holomorphically connected to the solution at p = with the condition c = Cm = 
/o + /i + 2 + 2m (m G Z>o). 

Now, we set 

exp(7rv^c™x) U'tl" 0{x + ij{p)) 



(5.13) ^Ti^,P) = A^(x,p) - i-iy^'Ami-x,p). 

By a similar argument to the case of the function A^"^{x), the function A'^"^{x,p) 
is square-integrable on the interval [0,1]. The function Af}("^{x,p) converges to the 
function Bl^Aff!^{x) (c = /o + ^i+2 + 2m) uniformly on compact sets for some non-zero 
constant B'^ as p ^ 0. 

Combined with the relation Aff!"^{x) = Bm^{x)ilJrn{x) for some B^, the following 
theorem is proved: 

Theorem 5.1. Let m be an non-negative integer. For each m, there exists e G 
]R>o such that if \p\ < e then there exists an eigenvalue Emip) and an eigenfunction 
gm{x,p) of the Hamiltonian of the Inozemtsev model (see (^. 3^) ) with the condition 
h = h = such that 

(5.14) Em{p) n^lo + /i + 2 + 2mf - — (/o(/o + 1) + /i(/i + 1)) and 

(5.15) gm{x,p) $(a;)^/'m(x) on compact sets, 
as p ^ and gm{x,p) is square-integrable, i.e. 

(5.16) / \gm{x, p)\'^dx < oo. 



Here, the function ipm{x) is defined in ( j5. 6| j. The functions Em{p) and gm{x,p) are 
holomorphic in p if \p\ < e. 

Remark For the case h = h = h = 0, this theorem was established in 0. □ 

5.3. The algorithm of the perturbation. Assume h = h = 0. From formula 
( |A.9| ), we obtain the following expansion in p, 

oo 

(5.17) h = Ht-Ct + J2 Mx)p'^ 

k=l 

where H is the Hamiltonian of the BCi (one-particle) Inozemtsev model (see (ST^)) 
with the condition l2 = h = 0, Ht is the Hamiltonian of the BCi Calogero-Sutherland 
model ( see ([4. 31) ) , Cj* is a constant, and Vk{x) are the functions which are expressed 
as a finite sum of cos27™2; in G Z>o). 
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Set Vm = then Vm is an eigenvector of the operator Ht with the 

eigenvalue Em = 7r^(/o + h + 2 + 2m)^. 

We now apply the algorithm of the perturbation. More precisely, we determine 
the eigenvalues Em{p) = Em + Yl^=i^m^P^ the eigenfunctions Vm{p) = Vm + 
YlkLi '^m' 4!^!m''^rn'P^ of the Operator Ht-Ct + YlkLi ^k{x)p'' as formal power series 
in p. 

Since the function Vk{x)vm is expressed as the finite sum of Vm' ijn' G Z>o) and the 
eigenvalues Em {m G Z>o) are non-degenerate, the coefficients of the formal power 
series Vm{p) and Em{p) satisfying 

(5.18) (^Ht - + £ M^)P'^ ^m{p) = Em{p)Vm{p), 

{Vm{p),Vm{p)) = {Vm,Vm), 

are determined recursively and uniquely (for details see |T0| ). 

The convergence of the formal power series Vmip) and Emip) is not guaranteed 
a priori and they may diverge in general. Nevertheless it is shown that the formal 



power series Vm{p) and Em{p) converge in our model as a corollary of Theorem ^TL 

Corollary 5.2. Let Vm{p) be the eigenvector and Em{p) the eigenvalue of the Hamil- 

tonian H of the BCi (one-particle) Inozemtsev model with the condition I2 = h = 0, 
which are obtained by the algorithm of the perturbation. If \p\ is sufficiently small then 
the power series Em{p) converges and the power series Vm{p) converges on compact 
sets of X. 

Proof. From Theorem p.l| , for each m there exists e G ]R>o such that if \p\ < e then 
there exists the eigenvalue Emip) and the square-integrable eigenfunction gm{x,p) 
that converge respectively to the trigonometric eigenvalue and eigenfunction asp ^ 0. 
Set gmix,p) = sj ^"^^^j^^irtiSy 9mix,p), then Em{p) and gm{,x,p) are holomor- 
phic in p near p = and they satisfy 

(5.19) - Ct + Yyk{x)p^ 9m{x,P) = Em{p)gm{x,p), 

{gmix,p),gmix,p)) = {^{x)i!mix),^{x)i)mix)) = {Vm,Vm)- 

These equations are the same as ( |5.18| ). From the uniqueness of the coefficients in p, 
we obtain Emip) = Em{p) and Vm{p) = gm{x,p). From the holomorphy of Em{p) and 
gm{x,p), the convergence of Emip) and Vmip) in P is shown. □ 

6. Comments 



6.1. After obtaining Theorem we should discuss the completeness of functions 
gmix, p) (see Theorem |5.1| ) and clarify properties of eigenvalues and eigenstates. These 
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matters will be discussed in the forthcoming article [|T0|. Note that some results will 
be obtained by applying the method performed in [Q. 



6.2. Roughly speaking, the problem of finding square-integrable eigenstates for the 
Hamiltonian of the Inozemtsev model (see ( f^.lD ) is equivalent to the problem of 
finding the solutions to the Heun equation (see (|2.4|)) with special conditions for a 
monodromy matrix around two singular points. 

It would be desirable if the Bethe Ansatz method would be of help in investigating 
the monodromies of the Heun equation. 
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Appendix A. 

This appendix presents the definitions and formulae for elliptic functions. 
Let uoi and uj^ be complex numbers such that the value uo^/^i is an element of the 
upper half plane. 

The Weierstrass p-function, the Weierstrass sigma-f unction and the Weierstrass 
zeta-function are defined as follows, 



(A.l) p{z) = p{z\2ui,2u-i) 



— + V 

z"^ ^ \{z — 2wiO\ — 2nu;3)^ {2'mu)\ + 2nijj-i]^ 

(m,n)eZxZ\(0,0) ^ V i i il 

o{z\ = p(z\2uJi, 2uj'i) = z TT I 1 I 



(m,n)6ZxZ\{{0,0)} 

■ exp 



C(^) 



2muJi + 2nuj3 2{2muji + 2nuj'iY J ' 
a'iz) 



a{z)- 
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Setting ijj^ = —uji — u-s and Cj = p{uji), rji = ({uji) {i = 1, 2, 3) yields the relations 

(A.2) ei + 62 + 63 = ?7i + ?72 + ?73 = 0, 

p{z) = -C'(x), {p'{z)Y = 4{p{z) - e,){p{z) - e,){p{z) - 63), 
p{z + 2iUj) = p{z), Ciz + 2LUj) = Ciz) + 2r]j, (j = 1,2,3), 



ip'{z)y 2\z-ei 2; -62 z-es 
p{z + uj,) = 6, + }^ ^, {% = 1, 2, 3 , 

where G {1,2,3} with i' < i", i 7^ i', and i 7^ i". The co-sigma functions ai{z) 
{i = 1,2,3) are defined by 

(A.3) ai{z) = exp{-r]iz)a{z + uji)/a{uji), 

then we obtain 

(A.4) (^^y = p(^)_e„ (z = 1,2,3). 

Let T be an complex number on the upper half plane. The elliptic theta function 
9{x) is defined by 



(A.5) 9{x) = 2 ^(-1)"-^ exp(r7rv^(n - (1/2))^) sin(2n - 1): 

n=l 

Then 

(A.6) e{x + 1) = -e{x), e{x + r) = -e-^^'^^-^^^eix). 

Setting T = LUs/uJi yields the relation 

(A.7) a{2ujix) = 2uJiexp{2r]iuJix'^)e{x)/e'{0). 



The following formulas are used to obtain Theorem 3.13: 



where ui = 1/2 and c<j3 = r/2. 

Setting c<Ji = 1/2, UJ3 = r/2, and p = exp(27r-\/^r), the expansion of the Weier- 
strass p function in p is given as follows: 

(A.9) p(^) = — 97 — T Stt^) — (cos2n7rx - 1). 

^ ^ ^ ^ ^ sin^ TTx 3 ^ 1 - p"^ ^ 

^ ' n=l 
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